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Abstract 



\o : 

0^ ■ 

I Short pulses of a probe laser have been used in the past to measure whether a two-level 

atom is in its ground or excited state. The probe pulse couples the ground state to a third, 
^ I auxiliary, level of the atom. Occurrence or absence of resonance fluorescence were taken to 

i—i • mean that the atom was found in its ground or excited state, respectively. In this paper 

I we investigate to what extent this procedure results in an effective measurement to which 

the projection postulate can be applied, at least approximately. We discuss in detail the 
complications arising from an additional time development of the two-level system proper 
during a probe pulse. We extend our previous results for weak probe pulses to the general case 
and show that one can model an ideal (projection-postulate) measurement much better with 

■ a strong than a weak probe pulse. In an application to the quantum Zeno effect we calculate 
I the slow-down of the atomic time development under n repeated probe pulse measurements 

■ and determine the corrections compared to the case of n ideal measurements. 
^ ■ 

0\ ; PACS numbers 42.50.-p; 32.80.-t; 03.65.Bz 

I . 1. Introduction 

C ■ 

■ 

^ I For a quantum system a measurement in general involves a macroscopic apparatus which 

O^' interacts vi^ith the system and leaves a pointer in some definite position. By an ideal measurement 

^ . at time t we mean a measurement whose effect on a state can be described by applying the 

^> I projection postulate at time t. The projection postulate as commonly used nowadays is due to 

\^ • Liiders |||]. For observables with degenerate eigenvalues his formulation differs from that of von 



Neumann The projection postulate has been widely regarded as a useful too^. It is also 
known that one can envisage more general measurements where the projection postulate is not 
applicable; cf. e.g. [^. 

Some time ago, Cook |Q] proposed to measure the ground or excited state (stable or metastable) 
of two-level a system by means of a short pulse of a probe laser which pumps the transition 
between the ground state, level 1, and an auxiliary rapidly decaying third level (see Fig. 1). 
After emission of a photon the atom will be in its ground state, and so it is natural to assume 
that occurrence of resonance fluorescence means that the atom is in level 1 at the end of the 
pulse. Absence of resonance fluorescence was assumed to imply that the atom is in level 2. This 
idea was subsequently used in an experiment of Itano et al. to test the so-called quantum 
Zeno effect (QZE) 1). 

^e-mail: beige@theorie.physik.uni-goettingen.de 
^e-mail: hegerf@theorie. physik. uni- goettingen.de 
^e-mail: sonderma@theorie.physik.uni-goettingen.de 

^Liiders stressed its provisional character: "The projection postulate will be employed only until a better 
understanding of the actual measurement process has been found" (G. Liiders, private communication to G.C.H.). 
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In a recent paper ||^] two of the present authors have investigated to what extent a short 
probe pulse as in Refs. [Q, ||] can be regarded as an effective measurement to which the projection 
postulate can be applied, at least approximately. We distinguished two cases. If there is no 
interaction between levels 1 and 2 during a probe pulse, we showed that a single probe pulse 
does indeed project on either |1) or |2), up to terms decreasing exponentially with its duration. 
If one has an additional small driving field (Rabi frequency O2) between levels 1 and 2 during 
a probe pulse the situation is more complicated, and we had to restrict ourselves to weak probe 
pulses (with Rabi frequency f^s much smaller than the Einstein coefficient A3 of level 3, but still 
^2 ^ ^3)- Under this additional assumption we were able to show that an atom with photon 
emission during the probe pulse is projected onto a state (density matrix) extremely close, but 
not identical to, and that an atom without photon emission is projected onto a state 

very close to |2)(2|. The small difference to the ideal projection-postulate result was explicitly 
calculated to first order in the small parameters 

We also pointed out that these e parameters had to be much smaller than 1 for a pulse of the 
probe laser to act as an effective measurement. For many probe pulses ("measurements") this 
small difference adds up. The cumulative effect was explicitly determined, and we used our 
results to analyze the experiment of Ref. |^] on the QZE. 

The QZE Q deals with rapidly repeated measurements on a system. Under the assumptions 
of n ideal measurements, at times At, . . . , nAt = t, t fixed, the QZE predicts a slow-down of 
the time development and ultimately a freezing of the state for At 0, or n 00. The QZE 
has found a tremendous interest in the literature 

In the experiment of Ref. Q several thousand ions are stored in a trap and an rf field 
in resonance is used to drive the transition between two levels 1 and 2, which are essentially 
stable. The rf field is a so-called vr pulse of duration T^^ which changes state |1) into state |2). 
The experiment then endeavors to perform n measurements of the populations of levels 1 and 
2 between and by means of n short pulses of a probe laser (see Fig. 2) which pumps the 
transition between level 1 and the auxiliary rapidly decaying level 3, with occurrence of resonant 
fluorescence from an atom taken to mean detection of the atom in level 1 and absence thereof 
as detection of the atom in level 2, as proposed in Ref. At time the actual population 
of levels 1 and 2 are then determined through the resonant fluorescence of all atoms in the 
trap under the influence of a long probe pulse, with the vr pulse switched off. The populations 
determined in this way are in good agreement with the predictions of the QZE . 

The relevance of this experiment for the QZE has been contested in the literature. It was 
pointed out by some authors [^, |ll], 12, 14, 15] that there is no need to use the QZE 



to explain the experiment. One could simply incorporate the auxiliary level and the driving by 
the probe laser together with the rf field in the Hamiltonian or in the Bloch equations of the 
corresponding three-level system. From a numerical solution of the Bloch equations one does 



indeed obtain the final population at time T^^ in good agreement with the experiment [11|. The 
probe pulses are, in this view, not measurements but part of the dynamics. 

In our previous paper we used our results on the close, but not perfect, resemblance of 
individual probe pulses with ideal measurements to determine in closed form, to first order in 
the small e parameters, the population of level 2 after n probe pulses, at the end of the vr pulse. 
This was compared with a numerical solution of the three-level Bloch equations, and an amazing 
agreement was found. In Ref. Q our conclusion was that the projection postulate is a useful 
tool to give a quick and approximate understanding of the type of experiment performed in 
Ref. , but that for a more precise description one needs a more detailed analysis of the actual 
measurements, as for example carried out by us in Ref. Q. 
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In the present paper we treat the general case of arbitrary, not necessarily weak, probe 
pulses. The duration of a probe pulse, Tp, is assumed to be less than = tt/Q2, and our 
results will be given to first order in the above parameters ep, and ea- It will turn out that 
with increasing 1^3, i.e. increasing strength, the probe pulse models an ideal measurement with 
projection postulate much better than for weak pulses; this is particularly true in the presence 
of the rf field. 

The plan of the paper is as follows. In Section 2 we deal with a single short probe pulse 
during which the rf field is also switched on, and the description of individual atoms is studied 
according to whether they do or do not emit photons during a probe pulse. We show that at 
the end of a probe pulse an atom can only be in one of two states, denoted^ by and p^, 
respectively, depending on whether or not the atom has emitted photons. These states still 
contain a contribution from level 3, which will decay during a short transient time. During this 
decay time the rf field will also change the state, and the state will depend on the particular 
decay time chosen. We show, however, that there are two uniquely defined states which effec- 
tively describe the resultant action at the end of a probe pulse and which one can regard as 
the projected states of an atom. The two states are independent of the initial state, have no 
contribution from level 3 and are independent of any particular choice of transient decay time. 
These uniquely defined states are, however, never exactly realized by the atom and are in this 
sense virtual and a mathematical tool. 

In Section 3 we consider repeated probe pulses in connection with the QZE and the cu- 
mulative effect of the corrections to the projection postulate. Closed expressions are given and 
compared to numerical solutions of the three-level Bloch equations. Our conclusion is that if one 
would perform the experiment of Ref. |p with a much stronger probe pulse (with Rabi frequency 

of the order of the Einstein coefficient A3) then the result for the level population would be 
much closer to the case of ideal projection-postulate measurements than for weak probe pulses, 
as in the original experiment of Ref. Q. 

In Section 4 we discuss our results, and in the Appendix we briefly review the quantum jump 



approach 17, 18, 19] which we use in this paper. The quantum jump approach is essentially 



equivalent to quantum trajectories [BO] and to the Monte-Carlo wave function approach 



2. Effect of a single probe pulse on the atomic state 



In this Section we discuss how a single probe pulse acts as an effective measurement of an 
atomic level and to what extent the usual projection postulate can be applied. 

The case where the rf field is switched off while the probe pulse is on was discussed in some 
detail in Ref. It was shown that if the length of the probe pulse Tp satisfies the condition 

Tp » max{A^\ As/nj} , (2) 

then the probe pulse provides an effective reduction of the initial state to |1) or |2) depending 
whether or not photons were emitted (allowing, in the former case, for a decay of the third- 
level contributions). The respective probabilities are the same as predicted by the projection 
postulate for an ideal measurement of the level populations. 

If the rf field is switched on during the probe pulse it causes a small pumping between level 
1 and 2. However, we will show that, independent of its state at the beginning of a probe pulse, 
right at the end of the pulse an atom will be in one of two states, fP and p^, depending on 
whether or not it has emitted photons. To allow for the decay of the third-level contributions a 
transient time at least of the order of several decay times A^^ is needed. During this time the rf 

tilde on a density matrix indicates that its trace is 1. We will also have occasion to use non-normalized 
density matrices. Here, a state can be pure or mixed. 
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field is also active, and therefore the result has a small dependence on the transient time chosen. 
At the end end of this section we therefore mathematically construct two density matrices, pp 
and pp, in the 1-2 subspace which are independent of the transient time. These states, although 
not physically realized, can be used in a consistent way as hypothetical or virtual states of an 
atom at the end of the probe pulse, i.e. as states on which the probe pulse projects. 



2.1 Subensemble without photon emission 



We now determine the state at the end of the probe pulse for an atom which has not emitted 
photons. The result is the same as that found in our Ref. 0], but the present derivation is more 
streamlined. According to Eq. (^) of the Appendix, an atom without photon emission evolves 

with [|i6|, m, m 

U,{^{t,0)=exp{-[H},^t/h} (3) 
where the reduced (or conditional) Hamiltonian is given by Eq. (^) of the Appendix, 

iH^^Jh = M = Mq + Mi (4) 



with 



Mo = il73{|3)(l| + ll)(3|}+ ^^3|3)(3| (5) 
Ml = ^172{|2)(1| + |1)(2|}. (6) 



The eigenvalues of Mq are = and 



lAs[l±^l-4nl/Ajj . (7) 



For M = Mq + Mi these go over into Aj, i = 1,2, 3, which are the roots of the cubic equation 

p{X) ^ det(M - A) = A^ - i^sA^ + ^{n^ + nl)x - ^A^nj = . (8) 

Newton's method or a linear expansion givesQ 

A2= ^02ep(l + 0(62)) (9) 

where O(e^) denotes a correction at least quadratic in e = (ep,e/j), so that 0(e^)/|ep remains 
bounded. For A13, one obtains 

X,,, = Xl,{l + 0{e)) . (10) 

If the Aj are different (which here means Ai 7^ A3) then M has three normalized eigenvectors 
|Aj), which are in general not orthogonal. The reciprocal basis vectors |A*) are defined by 

{Xi\X^) = Sij (11) 

and satisfy the eigenvalue equation 

Mt|A*) = Ai|A*) . (12) 

^In our case A2 is real and one has the exphcit bounds 

i J72ep(l + ei)-' ^ ^2 < i n2ep(l + e^)(l + 4)"' 
since p{X) changes sign between these bounds, as seen by insertion. 
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In particular, an elementary calculation gives 

IA2) = -iep|l) + i2)-e^j|3) + 0(e2) (13) 
|A2) = iepjl) + |2)-e^|3)+(!)(e2) (14) 

For e~*^* one has the representation^ 



3 3 
-Mt _ „-Ait I \ \ / \ j I _ „-Ait 



^e-^-*|A,)(A*| = ^e-^'*P, (15) 



i=l 1=1 



with Pi = |Aj)(A*|. If Tp satisfies the condition of Eq. (|^ then the exponentials exp{— Ai^aTp} 
can be neglected, while exp{— A2Tp} « 1. Therefore for any initial pure state \ip) one has 

e-*^^p|^) = e~^2^p(A^|?/;)|A2) + exponentially small terms (16) 

which is proportional to |A2). Thus at the end of a probe pulse the subensemble without photon 
emissions is, after normalization and up to exponentially small terms, in the state IA2) and thus 
independent of the initial state. 

The probability Po^Tp] lip)) for no emission during a probe pulse, for initial state is the 
norm-squared of Eq. ([l^). If one has a density matrix p = Oii\^pi) {ipil at the beginning of a 
probe pulse the no-emissions probability becomes 

Po{rp;p) = J2aiPo{Tp;\i>i)) . (17) 

i 

By Eqs. (|l^) and (|l4|) one obtains 



Po(rp;p) = P22-epvr^P22 + 2epImpi2-2e/jRe/)23 + C(e^) • (18) 



We thus arrive at 



p\Tp-p) ^ e-Mr, p^-M^r, ^ (Tpi />) | A2) (A2 | 

and, for the basis [1), [2), |3), we obtain from Eq. (p^) 

/ -iep \ 
pO = |A2)(A2| = iep 1 -eR + 0{e^) . (19) 
\ -en J 

We have thus shown that, up to exponentially small terms, a probe pulse projects each atom 
with no photon emission onto the state IA2), which is close to |2). The probability for this is 
given by Pq{tp; p) of Eq. (18) where p is the density matrix of the atom at the beginning of 
the probe pulse. In case the rf field is switched off {Q2 = 0) during the probe pulse, one has 
IA2) = |2) and the probe pulse projects onto |2) with probability P22, up to exponentially small 
terms. 



2.2 Subensemble with photon emissions 



If an atom does emit photons during a probe pulse, the emissions can occur at random 
times. Therefore, even if one had started with a pure state the subensemble of all atoms 

^Alternatively, 

e~^^* = — ^^c^^^* + cyclic permutations , 

(Ai — A2)(Ai — A3) 

which is directly checked by application to the eigenvectors. Comparison with Eq. ( p^ gives explicit expressions 
for I Ai) (A' I . In the above formula one can take the limit A3 — > Ai . This leads to a term with the factor t exp{— Aif}. 
The following arguments are then easily adapted to this degenerate case. 
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with emissions is described by a mixture. The corresponding density matrix is denoted by 
P^{t,p), with p the density matrix at the beginning of a probe pulse, and the normahzation 
of p-^ is chosen such that tTp'^{T,p) is the probabihty for emissions until time r. The aim of 
this subsection is to show that />^(Tp,p) is in general proportional to a p-independent matrix. 
Hence the subensemble with photon emissions is, at the end of a probe pulse, described by a 
fixed normalized density matrix p^ = p^/trp^, which is independent of the initial state. This 
density matrix p^ , which still contains contributions of level 3, will now be explicitly determined. 

We denote by I{t'; p) the probability density for the emission of a photon at time r'. Now, 
after a particular atom has emitted its last photon before r, at time r' say, the atom is in its 
ground state |1) and until r the time development is given by 

e-A/(--')|i)(i|e-A^n-r') = pO(^ _ (20) 
Hence the subensemble is described by 

P>{t-p) = [ dT'l{r'-p)p\r-T'-\l)l 

dT7(r-T';p)pO(T';Il)) . (21) 







This is also seen directly from Eq. (^) of the Appendix. The density matrix of the complete 
ensemble is given by 

p{t-p) = p\t-p) + p>{t-p), (22) 



and 



iTp>{T-p) = l-iTp\T-p) = l-Po{T-p) . (23) 



Since in our case /(r; p) = A^p^^{t\ p) by Eq. (|69|) of the Appendix, Eqs. (^) and (22) give the 
integral equation 

/(r; p) = A^pI^{t- p) + A, [ dr'lir - r'; p)p'i,{T'; |1)). (24) 

Jo 

This can be solved by Laplace transform]^. Since, by Eq. (pO|), p^^ is a sum of exponential terms 
whose exponents do not depend on p, I is of the form 

9 

I(t;/5) = co + cie-^i"+ E Cae"'^"" (25) 

a=2 

where Cj = Ci{p). One has Re/ij > 0, pi is real and of the order of A2, while all the other p^s 
are of the order of Ai and A3 0. 

Therefore one sees that cq is the stationary emission rate of the three-level system. Further- 
more, considering times of the order of the length of the probe pulse, i.e. exp{— /Uir} ~ 1 and 
Re //qT S> 1 for a > 2, one sees that cq -|- ci must be positive. Physically the emission rate for 
such times cannot exceed the stationary emission rate ^31^3/(^3 -|- 2^3) of the 1-3 system with 

= and initial state |1), and the same is true for very large r, for which exp{— ^ir} <C 1. 
Hence one has the inequalities 

0<( '^^I<-T#^• (26) 
~ 1 Co + Clip) - Al + 2nl ^ ^ 



^Alternatively one could solve the Bloch equations for p. 



^For $^2 = one has — 0, and in general /^i is quadratic in f22. The Hi's are in fact the eigenvalues of the 
matrix in the Bloch equations. Two of the eigenvalues vanish for ^2 = 0. From Eq. ( p^ one can determine /ii 
explicitly as ^ii = 2epn2(Ai + fll)/{Ai + 2!^i) + 0(^1). 



6 



Therefrom one obtains 



CO, \ciip)\ < ■ (27) 



Only these rough inequahties will be needed^. 
From Eq. (^) for e~*^* one obtains 



p°(r;|l))= ^P,|l)(l|P;e-(^'+^^)- (28) 



where Pi was defined as |Aj)(A*|. Inserting Eqs. (p5|) and ( p8| ) into Eq. (21) with r = Tp gives 

p>(rp;p) = (co+cie-'^^-p) H dr' p^r' ; 

Jo 

+ci jj^"' dr' (e-''^("p-"') - e-^i"p) p°(r'; [!)) 
+ T'dr' y c„e-'^"("p-"')e-2^2"'P2|l)(l|4 



a>2 



+ T'dr' y c„e-'^"("p-"') y e-(^'+^^)"'Pi|l)(l|P/ . (29) 

It will be shown that the third and the last term are of order e^. The last term is clearly 
proportional to exponentially small terms and can therefore be omitted. Since e~'^^^^p~'^ ^ — 
g-z^iTp < fj^^r' , the second term is bounded by 



cipi i ' dr' t'p\t';\1)) = cipi / dr' T'e"^"^'^^' P2\l){l\Pl 
Jo Jo 

+ci/ii J2 r'dr'T'e-(^>+^^)^'P,|l)(l|Pj . (30) 

(i,i)^(2,2) 

Now, P2II) = \X2){X^\1) is of order e, by Eq. (|l|), and P2|l)(l|-P2 is thus of orderQe^. The 
second integral in Eq. ( ^0|) is proportional to (Aj + Aj)~^. Multiplied by ci//i this becomes 
bounded by e^. Thus the second term in Eq. (^) is bounded by e^. For the third term the 
argument is again based on P2|l)(l|i^2 = C'(e^). Hence 

p>(Tp;p) = (co + cie-^i-p) r dr'p'ir'; |1)) + ©(e^) (31) 

Jo 

and thus one has for the normalized density matrix 

p>(Tp;/,)=p>(Tp) oc rdr' p\T';\l))+0{e') (32) 

JO 

which is, up to order e^, independent of the initial state p (cf. footnote 11). 



^"From Eq. (|2J) one obtains by a more detailed calculation 

1 AiQi 



Co 

= col -rT-^^-2-^-—^{X2\p\\2) 



2 -^4^ "4" ^3 

Normalization of involves division by trp^ = 1 — Po{Tp,p). For p close to |A2)(A2| this is of order e, and 
in this case Eq. holds only up to order tpT-^/rp. In the subsequent applications, however, this will play no 
role. 
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The integral can be computed in closed form by choosing for p in Eqs. and (^2|) the 
stationary state p^'' of the three- level system driven by ^3 and ^2- Then the emission rate 
/ = ^3^33 is time independent, and one obtains 

pss ^ ^-Mrp^.s^-Mtrp ^ ^^^.s ^0(^/. j^^) _ (33) 







With Eq. nsa) this yields 



p>(rp) = - e-^^p/^e-*^'^!') /tr(-) + ©(e^) . (34) 



One can express p** directly in terms of M a^^ 



(M-lA3)(Mt-i^3)/tr(.) . (35) 



One can understand Eq. ( p4D as follows. It takes much longer than the time Tp to reach the 
stationary three-level state, and therefore there is not enough time to build up an appreciable 
population of level 2. With the second term one just subtracts the excess population of this 
level from p^** because 

g-A/rp^sSg-A/trp ^ e'^^^rp ^^2 |^2^ |^^^ | ^ exponentially small terms . (36) 



From Eqs. (p4)'(p6D one can now calculate in a straightforward way the normalized density 
matrix p^{Tp) at the end of a probe pulse for the subensemble with photon emissions and obtains 



1 



= 42^002^ .2, o . ^iepO^Tp eR{Al + \ + 0{e^) . (37) 

^3 + 2i\ + A^epihTp \ ^^(^2 ^ j^2) j^2 

It is noteworthy that this is just, except for the e terms, the two- level stationary state of the 
1-3 system (with = 0). In addition p^ has a 22-component proportional to Tp. This results 
from the possibility of macroscopic dark periods for the V system under consideration [^] so 
that the emission of photons may stop before the probe pulse ends. Such atoms are then in a 
state close to |2) [|^] and contribute to P22- 

The effect of the probe pulse is therefore that right at its end an atom, which is initially 
described by />, is projected with probability 1 — Pq{tp; p) onto p^{Tp) and with probability 
Po{Tp; p) onto p^ . We note that p^ and p^ still have a 33 component which is small for pp but 
can be appreciable for p^ and whose decay will now be studied. 

2.3 Decay of the 33 component and the question of measurement by a probe pulse 

At the end of a probe pulse one still has 1^2 7^ while 0,^ = 0. The time development of an 
arbitrary density matrix is then governed by the corresponding Bloch equations with ^3 = 0, 



i^From Eqs. ^ and (|) one finds {M - ^A3)\3) = ^fisll) and M + Aft = ^3|3>(3|. From this one readily 
obtains 

dr e |l)(l|e 







"3^ 
4 



11^ A3 

Comparison with Eqs. (jj^ ) and gives ( ^sl) 



-(M4^3)f ..'A(e--e--')(M^4^. 
(m - 1^3) (1 - e-"-Pe-^^*-p) (Mt - 1^3) . 
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leading to a rapid decay of all i3 and 3i components, i = 1,2, 3. Simultaneously the 1-2 transition 
is weakly pumped by ^2- If P33 is very small the net effect of this is the same as increasing p-^^ 
by /jgg and replacing all pf^ and p^^ by right after the probe pulse. This was the case in Ref. 
1^]. If, however, p^g is not small, then the pumping from 1 to 2 during the decay is not quite 
the same with this replacement. 

Without employing Bloch equations one can very easily determine the time development by 
means of Eqs. ( |6^ and (|6^). Between probe pulses the reduced time development operator is 
given by exp{— Mfcr} with 



Mb 

Hence, with 

Un{T + To, To) 

and with the projector P12 = 



^1^2{|2)(1| + |1)(2|}+ ^^3|3)(3| . 



^ cos |r22T —i sin ^172''" 
= — isinir22T cos^^l2T 
^0 1 

|1)(1| + |2)(2| one has 
e-A^^- = C7,(r)Pi2 + e-^^^l3)(3| 



(38) 



(39) 



where describes the 1-2 pumping by the rf field and the last term the rapid decay of level 
3. Until time Tp + r the density matrix p^{Tp) has then developed, by Eqs. (|63D, ( |59D and ( |6^ ) 
of the Appendix, or by the analog of Eq. (^), to 



P>(Tp + T) 



e-^-^^-p>(rp)e-*<- 

+ Tdr' I{T';p>{T^))e 



'Mt{T-T') 



l)(l|e 



-MUt-t') 



(40) 



where, from Eq. (|3£ 

I{T';p>{T^)) = A,e-^^^'pUr^) 

An analogous equation holds for the subensemble fP{Tp + T) of Subsection 2.1. During the decay 
of the 33 component an additional photon may be emitted. Eq. ( ^ ) is easily evaluated for times 
r larger than a transient time rtr for which exp{ — ^^3Ttr} is at most of the order of and thus 
can be neglected. E.g. for 

Ttr = 15/As (41) 
the exponential is less than 10~^. Eq. (Hol) becomes, for r > rtr, 



P>(Tp + T) 



--12/0' (Tp)Pl2 



(42) 



With Eq. (^) for U-n- the integral is readily evaluated. Denoting the matrix in the brackets by 
Pp, one obtains 

^ PlliTp) + PssiTp) Pl2iTp) - ^P33iTp)(^A 
P21 l^p) + h P33 i^p)^A P22 i^p) 



p^ 



(43) 



V 















Insertion of p^{Tp) leads to 



P^ 



Al + inj + yl§epil2T, 



-iepAl + ^eA^l AlepQ.2Tp 




(44) 
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and thus, for r > Ttr, 

P>{t^ + t) = U^{t) p>uI{t) (45) 

up to the same order in e as /5^(Tp) (i.e. up to order e^, unless one had started close to |2) at 
the beginning of the probe pulse). 

In a similar way, one finds for the normalized density matrix of the subensemble without 
emissions during the probe pulse the density matrix at time Tp + r, for r > rtr, 

~p'ir, + T) = UAr) p'pUliT) (46) 

with 

/5p = 1^ iep 1 I , (47) 

again up to the same order in e as p^{Tp). 

These results can be viewed in the following way. At the end of a probe pulse both the 
subensemble with and without photon emissions still have components involving the third level, 
which subsequently decay during the transient time. Simultaneously the population of level 1 is 
increased correspondingly. In addition, during this decay, one also has the action of the rf field 
which introduces a change dependent on the length of the transient decay time one considers. 
However, from Eqs. (45) and (^) it is apparent that for times larger than the transient time 




one obtains the correct result also if one projects onto the states pp and pp at the end of a 
probe pulse and then develops with the rf field only. These projection states can be considered 
virtual in the sense that they are really never quite realized in the actual time development of 
the atom. But it is at least formally consistent to say that a probe pulse acts as an effective 
measurement and projects onto pp or pp. 

If one neglects all e terms the states pp and pp become |1)(1| and |2)(2|, respectively, as for 
an ideal measurement to which the projection postulate can be applied. The e terms can thus 
be viewed as corrections to the projection postulate. In the case <C 1 {0,2 <C ^43) the second 
term in the 1-2 component of pp can be neglected, and pp then reduces to Eq. (61) of Ref. 

3. Applications to the Quantum Zeno effect: Corrections to the projection-postulate 
results 

As shown above a probe pulse, regarded as a measurement pulse, models very closely an 
ideal measurement to which the projection postulate can be applied. As pointed out in Ref. [Q], 
for repeated probe pulses the small differences may add up, however, and the net result for an 
ensemble of atoms was calculated there in the case <^ A3. Here this will now be calculated 
for an ensemble of atoms without this restriction for n probe pulses during a vr pulse of length 
Tn = IT / 0,2 (cf. Fig. 2). The time between two pulses is denoted by AT. It will be assumed 
that AT ^ ^3 ^) in fact AT larger than the transient decay time rtr will do. 

At the end of the k-th probe pulse the total ensemble of atoms consists of two subensembles. 



one with and the other without emissions during the k-th pulse. Hence, by Eqs. ( |37D and (1£), 
the total density matrix p{t) at time t = k ■ (AT + Tp) is of the form 

p{t) = a{k)p> + f3{k)f + 0{e^) (48) 

with a + (3 = 1. We will now determine a{k) and f3{h). 

We consider an ensemble with density matrix at the end of a probe pulse, develop for the 
time AT without the probe pulse and then again switch on the next probe pulse. We define p as 
the probability of finding no photons during this pulse. Similarly, q is defined as the no-photon 
probability for starting with fP instead of fP" . We will calculate p and q explicitly further below. 
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Now, if the density matrix is given by Eq. ( [48D at the end of the k-th probe pulse, then at 
the end of the next pulse it is given by 

(1 - p)a{k)p> + paik)p° + (1 - q)P{k)p> + qPik)p^ . (49) 

Hence f3 satisfies the recursion relation 

P{k + 1) = paik) + qp{k) . (50) 

Using a = 1 — (3 the solution is seen to be 

for k > 1. We determine /3(1) for the case that at t = all atoms are prepared in the ground 
state. At the beginning of the first probe pulse the state is then C/^(AT)|1), and /3(1) is the 
probability for no photon emission until the end of the first pulse. With Eqs. ( |3^ ) and (|l^) and 
the abbreviations 

c = c(AT) = cos(J)2Ar) , s = s{AT) = sin(02Ar) (52) 

one easily obtains 

Pil) = 1(1 - c) + eps - ^ vr ^ (1 - c)ep + 0{e^) . (53) 

In the determination of the no-photon probabilities p and q the relevant density matrices at 
the beginning of the next probe probe pulse are [/^(AT) pp U.n-{AT)^ and UTr{AT) pp Ut^{AT)^ . 
Using Eq. (18) one obtains, with eA = 1^2/^3, 

q = 1(1 + c) - 6p |2s + 1 f)2Tp(l + c)| + 0{e^) . (54) 

One can now insert /3(fc) of Eq. (^) and a{k) = 1 — /3(/c) with k = n into Eq. (^) to obtain 
the density matrix p{t = T^^) at the end of the vr pulse, after n probe pulses. Expanding in terms 
of e one obtains for the population of level 2 

P22(T.) - -(l-c)+6pj.c ^2 + 2172 

Ai + Oi l-c" / , rp\ 1 



+7r -f nc" -.^^ s + vr -f 



where c(Ar) and s(Ar) are given by Eq. (||), ep = ^302/^^3 < 1 and = ^22/^3 < 1- 

This "quantum jump" result contains additional e terms when compared with the special 
case in Eq. (76) of Ref. for <C A3, while up to zeroth order it is the same. The zeroth 
order gives the approximation 

MT.)^i{l-cos"(i-i)} (56) 
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n 


Projection Postulate 

AT = T-j^/n AT = T-j^/n — Tp 


Quantum Jump 


Bloch eq. 


Observed Q 


1 


1.00000 


0.99978 


0.99978 


0.99978 


0.995 


2 


0.50000 


0.49957 


0.49960 


0.49960 


0.500 


4 


0.37500 


0.35985 


0.36062 


0.36056 


0.335 


8 


0.23460 


0.20857 


0.20998 


0.20993 


0.194 


16 


0.13343 


0.10029 


0.10215 


0.10212 


0.103 


32 


0.07156 


0.03642 


0.03841 


0.03840 


0.013 


64 


0.00371 


0.00613 


0.00789 


0.00789 


-0.006 



Table 1: Predicted and observed population of level 2 at the end of the vr pulse for n 
probe pulses of length Tp . 





Projection Postulate 






n 


AT = T^/n - Tp 


Quantum Jump 


Bloch eq. 


1 


0.99978 


0.99978 


0.99978 


2 


0.49957 


0.49957 


0.49956 


4 


0.35985 


0.35985 


0.35979 


8 


0.20857 


0.20858 


0.20853 


16 


0.10029 


0.10030 


0.10027 


32 


0.03642 


0.03642 


0.03641 


64 


0.00613 


0.00613 


0.00613 



Table 2: Predicted population of level 2 at the end of the tt pulse for the parameters of 
the experiment, but with Q3 = ^3/2. 



while for n instantaneous ideal measurements the projection postulate would yield 

- (l-cos" -1 . (57) 



2 I n 



Thus the zeroth order result in Eq. ( [Sq ) can be viewed as the result of the projection postulate 
with the finite duration of the measurement pulse taken into account. The intuitive reason for 
this has been discussed in some detail in Ref. 0. 



In the following tables we have evaluated the expressions for p22{Tn) of Eq. (56), i.e. the 
(modified) projection-postulate result, and of Eq. (^5|), the quantum jump result, as well as a 
numerical solution of the corresponding Bloch equations of Eq. ( [70|) of the Appendix. 

Table 1 is calculated for the parameters of Ref. |], i.e. T^ = 0.256s, A3 = 1.2 • lO^/s, 
^3 = 1.9 • 10^/s, and Tp = 2.4 • 10~^s; here <C A^. The column with the quantum-jump 
results coincides with the corresponding column in Table 2 of Ref. Q. This shows that the 
additional e terms in Eq. (^) compared to Eq. (76) of Ref. are indeed negligible in this 
case. The agreement between the quantum-jump results and Bloch equations is excellent, while 
the column with the (modified) projection-postulate results show the small, but still markedly 
noticeable, relevance of the e terms in the quantum-jump expression of Eq. (]55|). In this case the 
probe or measurement pulses give a very good, but not perfect, realization of ideal measurements. 

A much better realization of ideal measurements can be obtained by choosing 0,^ much larger, 
e.g. of the order of A3. Then one has a strong pumping between levels 1 and 3, with many 
photon emissions, and the last photon during a probe pulse will therefore be emitted shortly 
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before its end. After this the atom is in the ground state, and the rf field has httle time to build 
up a correlation between levels 1 and 2. On the other hand, the population of level 3 at the end 
of a probe pulse grows with Og, and during the transient decay time after the end of a probe 
this leads to a build-up of the 1-2 correlation, but with the opposite sign, as seen from pp in 
Eq. (0). The net result is that the probe pulse projects the subensemble with photon emissions 
onto a state much closer to |1)(1| than for small 0- ^'^^ large ^3 one therefore expects the 
projection-postulate results to be much closer to the quantum jump or Bloch equations results. 
This is indeed borne out in Table 2 for 0,^ = A3/2, the other parameters being as in Table 1. 
Now the projection postulate results (with finite pulse duration) in the first column agree with 
the other ones much better, with the difference starting in the fifth decimal. 



4. Discussion 



In this paper we have investigated in detail to what extent a short pulse of a probe laser, 
which pumps the ground state of an atom to a third auxiliary level, can be regarded as a 
measurement of the states of the two-level system. In contrast to our previous paper |^ we have 
now allowed an arbitrary strength of the probe pulse. We have shown that in this general case 
the projection-postulate result is modified by additional small correction terms. But we have 
also shown that for a strong probe pulse, e.g. $^3 in the order of A^, the probe pulse acts in a 
way much closer to the ideal projection postulate. This is evident from the "virtual" projection 

which for large are extremely close to |1)(1| and 
1 2) (2 1, respectively. Indeed, the off-diagonal elements, which are already small for weak probe 
pulses, become orders of magnitude smaller for strong probe pulses. 

The slow-down of the time development of a state under repeated probe pulses ("measure- 
ments") is apparent in Tables 1 and 2. For strong probe pulses this is much closer to the 
projection-postulate predictions (with finite pulse duration taken into account). But an actual 
freezing of the state, as predicted by the projection postulate for instantaneous ideal measure- 
ment for At — > 0, is of course still not obtainable by these probe pulses. If one decreases the 
time AT between the pulse to the order of the transient time rtr or even to ^3 ^, then the third 
level does not decay completely to level 1, and the probe pulse can no longer be regarded as a 
measurement pulse. 



matrices pp and pp of Eqs. (|45|) and (46 



Appendix: The quantum jump approach in quantum optics 



We briefly summarize the quantum jump approach used in this paper. The quantum jump 
approach ||l^, 17, 18, 1^], quantum trajectories |2^ and the Monte-Carlo wave function approach 
iPUl are essentially equivalent. It describes a radiating atom between photon detections by a re- 
duced (or conditional) time evolution operator giving the time development under the condition 
that no photon has been detected [16|. After a photon detection one has to reset the atom to 
the reset state ("jump"), with ensuing reduced time development, and so on. The general reset 
states have been determined in Ref. [^]; cf. also Ref. [p!9|]. For a driven system with many 
emissions one then obtains a stochastic path, also called a quantum trajectory [20|. For a V 
system as considered in this paper the reset state after an emission is the ground state. The 
reduced time development together with the reset states provide a complete stochastic descrip- 
tion of the time development of the atom |18, 19 1. Starting with this description one can then 
derive the Bloch equations describing an ensemble of radiating atoms |16, |l^ ]. 

We consider the V system depicted in Fig. 1. In this system the upper levels 2 and 3 couple 
to a common ground level 1, with Einstein coefficient ^3 (in this paper level 2 is taken as 

^^Note that ep — Q,2As/fl^ decreases quadratically with Qs. Hence the off-diagonal elements in pp and pp 
become extremely small when is increased to the order of A3. 
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stable). We assume here that UJ32 = W3 — 1^2 is in the optical range. For simplicity we consider 
zero detunings of the driving fields, whose (real) Rabi frequencies are denoted by and ^3, 
respectively. In the interaction picture with respect to the free atomic Hamiltonian H^^, the 
reduced Hamiltonian H^^^ is given by ||^, 16, ^] 





hLJ^ =- 1^2 \ ^-iM. (58) 



where the atomic operator M is defined by the l.h.s. and where we have used matrix notation 
with respect to the atomic basis |1), |2), |3). The time development of an atom between emissions 
is then given by 

C//^^(r) = e-^^-d-/?i = e-^^" . (59) 
The no-photon probability until time r is then, for initial state 

Po(r;|V')) = ||e-^nV)f (60) 
or, more generally for an initial density matrix p, 

Po(r;/5)=tr{e-^V-"''"} • (61) 

The probability for the first photon to be emitted in (r, r + dr) equals Po(''"; p) — Po{t + dr; p) = 
wi{t] p)dT, where 

wi{t;\p) = -^ Po{t;p) (62) 

is the probability density for the first photonf^. For small upper level separation nonzero off- 
diagonal generalized damping terms may appear which lead to interesting coherence effects 



p3| , 24, 25, |26[]. For general n- level systems the reduced Hamiltonian is derived in Ref. |18| , |lSf| . 

The reduced time development is not unitary. The reason is that it does not describe the 
time evolution of the whole ensemble but that of the subensemble with no photons. The size 
of this subensemble is decreasing in time since an atom for which a photon has been detected 
leaves the subensemble, and this is refiected in the decrease of the norm squared in Eq. (|60|). 
The above probability density determines the (random) time for the first photon. After that the 
atom is reset, for a V system to the ground state, |1). The next emission time is then determined 
by wi{t; |1)), and so on. In this way one obtains a quantum trajectory. 

From this description of single systems one can recover the usual Bloch equations of the 
complete ensemble as follows The density matrix p{t) of the ensemble is a sum of two 

terms, p-^ and p^, corresponding to a subensemble of atoms with and without photon emissions 
until time r, respectively, 

^(r) = pO(r)+p>(r) . (63) 

From Eq. (^) one has for initial state p 

p°(r;p)=e-^^>e-*^^-. (64) 

If I{t'; p)dT' denotes the (unconditioned) probability to find a photon in (r',r' + dr') then, at 
time r, the sub-subensemble of atoms with their last emission in this interval is described by 

I(r';p)dry(r-T';|l)) (65) 



^■^Depending on the parameters there may be a finite probabihty that no photon is emitted at all. Therefore 
this probability density need not be normalized to 1. 
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and therefore 

P>(t;/9)= r dr' I{t';p)p\t-t'-\1)) . (66) 

JO 

Thus one has 

p{t) = e-^- p e-^'^ + r dr' I{t';p)p^{t - r'; |1)) . (67) 

JO 

Differentiation gives 

p(r;p)=pO(r;p)+/(r;p)|l)(l|+ f dr' I (r' ; p) p'' (r - r' ; . (68) 

JO 

Taking the trace and using tr p(r) = 1 gives 

I{T;p) = A3P33{r;p) (69) 
and thus Eq. ( p7| ) becomes in the present situation 

p(r;p) = e-*^>e-*^'^ + dr' A^pssir'; p)p'{T - r'; |1)) . (70) 

JO 

From Eq. (|m|) one obtains p'^, and inserting this into Eq. (|6^ ) gives 



/^(r) = -^[F,Xr) -p(t)f/J"J +A3P33(t)|1)(1| . (71) 

This is a compact form of the Bloch equations used in Refs. pO| , pT| ]. Conversely, from Eq. ( [7l| ) 
one can immediately obtain the integral equation of Eq. (^) . 
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Figure 1: V system with (meta-)stable level 2 and auxiliary level 3 with Einstein coefficient 
As. 0,2 and are the Rabi frequencies of the rf field and the probe laser, respectively. 



AT 



vr pulse 



probe pulse 



Figure 2: Probe pulses and n pulse 



18 



